Microwave-induced flow of vortices in long Josephson junctions 
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We report experimental and numerical study of microwave-induced flow of vortices in long Joseph- 
son junctions at zero dc magnetic field. Our intriguing observation is that applying an ac-bias of a 
small frequency / -C f p and sufficiently large amplitude changes the current-voltage characteristics 
(I-V curve) of the junction in a way similar to the effect of dc magnetic field, well known as the 
flux-flow behavior. The characteristic voltage V of this low voltage branch increases with the power 
P of microwave radiation as V 3 cx P a with the index a ~ 0.5. Experiments using a low-temperature 
laser scanning microscope unambiguously indicate the motion of Josephson vortices driven by mi- 
crowaves. Numerical simulations agree with the experimental data and show strongly irregular 
vortex motion. We explain our results by exploiting an analogy between the microwave-induced 
vortex flow in long Josephson junctions and incoherent multi-photon absorption in small Josephson 
junctions in the presence of large thermal fluctuations. In the case of long Josephson junctions 
the spatially-temporal chaos in the vortex motion mimics the thermal fluctuations. In accordance 
with this analogy, a control of the intensity of chaos in a long junction by changing its damping 
constant leads to a pronounced change in the shape of the I-V curve. Our results provide a possible 
explanation to previously measured but not yet understood microwave-driven properties of intrinsic 
Josephson junctions in high-temperature superconductors. 

PACS numbers: 74.50.+r; 74.25.-Nf; 85.25.Cp 



I. INTRODUCTION 

The discovery of intrinsic Josephson junctions in highly 
anisotropic high-T c superconductors [lj has boosted the 
interest in the dynamics of Josephson vortices. These 
junctions were proposed as candidates for high power 
flux-flow oscillators at THz frequencies. The Joseph- 
son penetration depth Aj in these junctions is rather 
small, typically of the order of 0.1 /im. Thus the stan- 
dard micron-size intrinsic junctions possess motion of 
Josephson vortices, similarly to conventional low-T c long 
(I ^> Xj) junctions. 

Many of the experiments with intrinsic junctions have 
been performed during the last decade at microwave 
(GHz) frequencies, which are well below the junction's 
zero-bias plasma frequency / p . These experiments led to 
observations of Josephson plasma resonance 2] and novel 
microwave-induced branches in the current- voltage char- 
acteristics (I-V curves) of intrinsic junctions 0. While 
the Josephson plasma resonance has been thoroughly in- 
vestigated and explained, the effects produced by mi- 
crowaves on I-V curves have not yet been fully under- 
stood. One obvious reason for that is that the dynam- 
ics of Josephson junctions at frequencies / < / p can be 
chaotic and thus complicated. Nonetheless, many exper- 
iments with underdamped intrinsic junctions driven by 
microwaves at zero dc magnetic field revealed very puz- 
zling smooth branches in I-V curves 3]- The aim of our 
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work is to find an explanation for this behavior. For 
this purpose, we have performed a detailed study of con- 
ventional low-T c long Josephson junctions under strong 
microwave irradiation and have done also numerical sim- 
ulations, which we compare with experimental data. 

The dynamics of long Josephson junctions under mi- 
crowave irradiation was studied in various contexts in the 
past both theoretically and experimentally. In particular, 
there had been an extensive search for regimes in vortex 
motion synchronized by the external microwaves . The 
vortex dynamics in underdamped junctions under the in- 
fluence of ac drive is often very complicated, showing 
along with synchronized (phase-locked) dynamics |5| also 
the chaotic behavior The spatially-temporal chaos 
in underdamped long Josephson junctions typically oc- 
curs at low drive frequencies (/ < / p ). Moreover, the 
chaotic behavior in long junctions may appear even in 
the absence of external irradiation at low characteristic 
frequencies of vortex oscillations |7| . Increase of damping 
usually helps obtaining more regular dynamics. 

The interest to drive vortices in long junctions by 
microwaves has been also biased by practical applica- 
tions. Recently, low-frequency microwave properties of 
long Josephson junctions have been studied in relation 
to their possible application in tuneable resonators 
On the other hand, driving long junctions by microwaves 
is now successfully used for phase-locking of Josephson 
flux-flow oscillators in integrated superconducting sub- 
millimeter wave receivers |9| . 

In this paper, we investigate experimentally (Section 
II) and numerically (Section III) the effect of microwave 
irradiation on long underdamped Josephson junctions. 
In particular, we are interested in the influence of mi- 
crowaves on the current-voltage characteristics of the 
junctions at zero (or small) dc magnetic field. Our 
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results suggest an explanation for previously observed 
microwave-induced smooth branches in I-V characteris- 
tics of intrinsic high-T c junctions. We find that these 
smooth branches are due to flux flow, but the vortex 
motion is very irregular and shows spatially-temporal 
chaotic behavior. 

It is interesting to note that the smooth branches in 
the I-V curves of long junctions resemble the ones ob- 
served in extremely small Josephson junctions driven by 
microwaves of a small frequency in the presence of large 
thermal fluctuations ^(j- I* 1 this case the effect of mi- 
crowaves on the I-V curves has been explained as inco- 
herent multi-photon absorption in the Josephson phase 
diffusion regime. Thus, we will discuss (Section IV) our 
results by exploiting the analogy between the spatially- 
temporal chaos of vortex motion in long junctions and 
the diffusive behavior of the Josephson induced by ther- 
mal fluctuations in small Josephson junctions. 
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II. EXPERIMENTS 



We performed measurements using Nb/Al-AlO^/Nb 
long Josephson junctions 1 1 1 1 ] of overlap geometry. Be- 
low we present our results for the junction of the length 
I = 250 /im and width w — 3.5 /im. The critical cur- 
rent density j c is about 180 A/cm 2 , which corresponds to 
the Josephson penetration length Aj of about 28 /im and 
zero-bias plasma frequency of f p ~ 77 GHz. The sam- 
ple was placed inside the vacuum volume of a variable- 
temperature optical cryostat. Magnetic held was applied 
in the plane of the junction perpendicular to its larger 
dimension. The dc measurements were carried out using 
a four-probe configuration. 

We investigated the junction behavior under the mi- 
crowave irradiation of a frequency in the range from 
10 GHz to 20 GHz at two different temperatures 4.2 K 
and 5.0 K. The external microwave radiation was applied 
to the junction by an antenna placed outside the cryostat 
close to the optical window facing the sample. Thus, the 
ac bias current was induced by the antenna in the bias 
leads and across the junction. The typical power P of 
the applied microwave radiation was ranging from few 
pW to several mW, as referenced to the output port of 
the microwave source. Below we present and discuss the 
results obtained for the frequency of microwave radia- 
tion / = w/(2tt) = 16.61 GHz, which corresponds to the 
/ ~ 0.2/p. We note that results obtained at different 
frequencies within the studied range were qualitatively 
similar. 

The dependence of the critical current of the junc- 
tion on microwave power at zero magnetic held and tem- 
perature 4.2 K is shown in Fig. da). The critical cur- 
rent I c monotonically decreases with power and drops 
abruptly to zero at the power level P cr ss 3 dBm. This 
kind of behavior can be interpreted in general terms 
for nonlinear systems as a threshold in energy trans- 
mission [T^ . For a long Josephson junction, it can be 
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FIG. 1: (a) Critical current I c of the junction at zero magnetic 
field versus the applied microwave power P at the frequency 
/ = 16.61 GHz. (b) I-V curves for different power levels at 
P > P cr . The measured voltage V a at the top of the first 
curve is also shown, (c) Dependence of the voltage V s on the 
microwave power P using double-logarithmic scale. 



3 



treated as ac-driven dissociation of a virtual breather 
(vortex-antivortex bound state) at the junction bound- 
ary A smoother drop occurs at the power level 
of P r w — lOdBm. It is associated with the Josephson 
plasma resonance, which we will not address in the rest of 
the paper. Plasma resonance in long junctions has been 
studied previously in several experiments Il5| . 

As the applied microwave power exceeds P cr , the su- 
percurrent branch of the I-V curve disappears. Fig.^b) 
shows the gradual changes in I-V curves with increasing 
microwave power above P cr . Each I-V curve has a linear 
region at low bias and its slope decreases with power. At 
a certain bias point, which depends on the power level, 
we observe a sudden jump from the voltage V = V s to 
much higher voltages of the energy gap at V « 2.7 mV. 
The dependence of the voltage V s of this critical point 
on the applied power is presented in Fig. |TJc). We have 
found this dependence to be power-law V s cx P a with the 
index a — 0.55 ± 0.01. Very similar behavior of the I-V 
curves was also observed at the temperature 5K. The 
major difference was in the modified shape of the I-V 
curves, where the linear branch at low bias was replaced 
by a more pronounced pseudo-resonant feature, i.e. a 
peak with downward curvature. Figure GJa) illustrates 
this behavior. The microwave power dependence of the 
switching voltage V s (its position is marked in Fig. Ufa)) 
remained power-law with a = 0.62 ± 0.01. 

It is worth mentioning that there is some uncertainty 
in determining the position of the switching voltage V s 
(especially in the case of numerically calculated curves - 
see results below). At the same time, it will be useful 
in the following to have a criterion for the voltage which 
characterizes the effect of microwaves on the I-V curves. 
One could, for example, choose some current value and 
observe how the voltage Vj at this current changes with 
microwave power P. In such a case, however, question 
arises whether the behavior of Vj(P) is universal or de- 
pends on current. We have investigated Vi{P) depen- 
dence for three current values marked by dotted horizon- 
tal lines in Fig.^a). The dependence is always power-law 
Vi oc P 13 at j3 = 0.68, 0.63 and 0.47 for current values 
150; 184 and 223 /iA correspondingly. Thus, the power- 
law V tx P x , where V is either V s or Vi and value of 
x is close to 0.5, describes the IV-curve evolution versus 
applied microwave power with a reasonable accuracy. 

It is interesting that the I-V curves under the influence 
of microwave irradiation resemble the well known flux- 
flow state of long Josephson junctions. The microwave 
field y/P affects the junction similarly to dc magnetic 
field H applied in the plane of the tunnel barrier. How- 
ever, the sudden drop of the junction critical current I c 
to zero at P = P cr qualitatively differs from behavior of 
I c in dc magnetic field. 

In order to experimentally obtain information about 
the microwave-induced vortex motion in the junction we 
have performed spatially-resolved measurements using a 
low-temperature Laser Scanning Microscope (LSM). 

The principle of probing the junction by laser beam is 




FIG. 2: Schematic view of the laser scanning experiment. 
Dimensions are not to scale. 



to introduce a local heating of the sample at the beam 
position and to get the information about the local prop- 
erties of the junction at this point. We scanned a focused 
laser beam over the junction area and simultaneously 
recorded its electrical response. The voltage response of 
a current-biased junction is due to the change of damping 
induced by local temperature increase at the illuminated 
spot. In order to improve the signal-to noise ratio and in- 
crease spatial resolution, the beam intensity is modulated 
at a frequency /mod and the circuit response is detected 
using lock-in technique. The spatial resolution of this 
technique is limited by both the thermal healing length 
(depending on the material properties and / mo d) and the 
laser beam diameter. In the presented experiment it was 
about 2 /zm. 

Figurc[2]illustrates the scheme of the LSM experiment. 
The junction area is scanned by tightly-focused laser 
beam, which power is modulated at / mo d = 100 kHz. 
The voltage response SV arising due to the tempera- 
ture oscillations is detected by lock-in technique. The 
characteristic time of a laser beam displacement is much 
longer than l// mo dj which allows to measure the volt- 
age response for every geometrical point of the junction 
averaged over many thermal cycles. In its turn, the char- 
acteristic modulation time l// mo d is much longer than 
the flux-flow oscillation frequency Vj $o in the junction, 
where $o = h/(2e) is the magnetic flux quantum. We 
adjusted the laser beam power to have the relative mag- 
nitude of the voltage signal |AV/V| — 10 -3 , which en- 
sured that the beam acts as a weak perturbation of the 
dynamical state of the junction. In order to reduce noise, 
we averaged the measured SV over the junction width. 
We plot the inverted voltage response (—SV) as a func- 
tion of the laser beam position along the length of the 
junction. 

Results of the LSM investigation of the junction state 
corresponding to point A marked on the I-V curve in 
Fig. HJa) are presented in Fig. OJb). The vertical axis 
displays the inverted voltage signal (—SV) in arbitrary 
units. Solid circles illustrate the LSM response in the 
absence of dc external magnetic field. There is hardly 
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FIG. 3: (a) I-V curves for different power levels obtained at 
the same regime as ones from Fig. 0^b) but at higher tem- 
perature T = 5 K; (b) LSM response of the junction in a 
state corresponding to the marked bias point A in (a), at zero 
(solid circles) and non-zero dc magnetic fields (H = 2.3 Oe up 
triangles, H — —2.3 Oe down triangles). 



any voltage response at the junction edges, while a no- 
ticeable maximum occurs in the center of the junction. 
The pattern is rather symmetric reflecting the symmetry 
of the junction at zero magnetic field. 

To check the above hypothesis of the flux-flow nature of 
the observed microwave driven behavior of the junction, 
we also measured LSM patterns at non-zero dc magnetic 
field directed perpendicular to both current flow direction 
and junction length as shown in Fig. [21 The main idea 
of applying dc field is to break the symmetry between 



two junction edges and facilitate injection of Josephson 
vortices at one of them. We would expect that under 
such conditions the vortices would be mainly entering 
the junction from one edge. The measured response pat- 
terns corresponding to two magnetic fields of 2.3 Oe equal 
in magnitude but opposite in polarity are presented in 
Fig. |HIb) by up- and down-triangles. We clearly observe 
asymmetric patterns: For H > a maximum of the re- 
sponse is located at the right half of the junction while in 
the case of H < the maximum is shifted to the opposite 
junction edge. We obtained similar asymmetric patterns 
in non-zero magnetic field for all measured bias points of 
I-V curves in the Fig. |3fa), as well as for bias points in 
Fig. ^b). The maximum in the voltage response —SV 
was shifting back and forth between junction edges, de- 
pending on the polarity of the applied field. This fact 
indicates that the response is most probably generated 
by Josephson vortices entering the junction from one of 
the edges. 

We note that similar asymmetric response patterns 
were observed for a pure flux-flow state in a long Joseph- 
son junction by Quenter et al. 0] and analyzed in 
Ref. 01- The response is typically the largest in the 
region where vortices are injected in the junction. In de- 
tail, the effect of a local perturbation induced by a hot 
spot (e.g., laser beam) on vortex dynamics in the junction 
has been studied elsewhere pcj . 

Thus, we argue that junction dynamics under the mi- 
crowave irradiation at relatively low frequency / -C fp 
is governed by the injection of Josephson vortices at the 
junction edges. This hypothesis is further supported by 
numerical simulations presented in the following section. 



III. NUMERICAL RESULTS 

The ac-driven weakly damped long Josephson junction 
is described by the well-known perturbed sine-Gordon 
equation for the Josephson phase difference ip(x,t) 



tptt - sin (p = a(x)ip t - 7 - 7 ac sin(wf ) , (1) 



where the subscripts stand for the partial derivatives, 
the coordinate x, running along the junction, and the 
time t are normalized so that the Josephson penetra- 
tion length Aj and plasma frequency u p — 2nf p are both 
equal to unity. A parameter a{x) is the damping coef- 
ficient, which we will assume to be spatially dependent 
when simulating the effect of the laser beam on the junc- 
tion. The term 7 is the uniform bias current density 
normalized to the critical current density j c , and 7 ac is 
the amplitude of uniformly-flowing ac current normalized 
in the same way. 

If the long junction is placed in external dc and ac 
magnetic field, Eq. JTJ) is accompanied by the boundary 
conditions: 



tp x (0,t) = h - /i ac sin(u;t), 



(2) 
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(fi x (£,t) = h + h ac sm(u;t), 



(3) 



where i = L/Aj is the normalized length of the junc- 
tion. In (J2J and © we take opposite signs in front of h ac 
assuming that a part of ac bias (due to its nonuniform 
distribution) flows near the junction edges [3 El- The 
precise distribution of the oscillatory currents and fields 
induced by microwaves in the junction can be rather com- 
plex. We suppose that in the absence of a ground plane 
the largest microwave current is induced near the junc- 
tion edges. That can be described by a distribution of 
the current 7 ac sharply peaked near the junction bound- 
aries or, equivalently, by an oscillating bias component 
± h ac at the edges. In the following we show the simu- 
lation results obtained for 7 ac = 0, so that microwave 
drive is coupled via the boundary conditions. We have 
also done simulations with a uniformly distributed 7 ac 
current and got qualitatively similar results. 

We numerically integrated Eq. Q with boundary con- 
ditions (0 and J3J|. The details of the simulation pro- 
cedure can be found elsewhere |20|. Simulations were 
performed with parameters very close to our experiment 
described above: a normalized junction length t = 8.9, 
the normalized ac-drive frequency u>/ (2ir) — 0.2, and two 
values of dissipation coefficient ot(x) = ao = 0.05 and 
a = 0.02. The calculated averaged voltage 8 is just 
(tp) normalized to the voltage at the first zero-field step 
(ZFS1) (^)zfsi = 2n/L 

Figure EJa,b) shows the numerically calculated 
current-voltage characteristic for several values of the 
normalized ac field amplitude h ac at zero dc magnetic 
field h = 0. As the amplitude of ac drive h ac approaches 
the amplitude h ac — 2, the critical current is reduced and 
the I-V curves display a smooth branch at low voltages. 
In this regime, in order to obtain accurate enough curves 
we had to perform numerical integration over 40,000 time 
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FIG. 4: Numerically calculated I-V curves for different am- 
plitudes of the ac bias and two different values of dissipation: 
(a) ao = 0.05; (b) ao = 0.02. Dashed IV-curve shows the first 
zero-field step (ZFS1) corresponding to a vortex shuttling be- 
tween junction boundaries in the absence of microwaves, (c) 
Dependence of the voltage Qj at the dc bias levels 7 = 0.072 
(circles) and 7 = 0.088 (triangles) indicated by horizontal 
dotted lines in (b) on ac-drive amplitude /i ac using double- 
logarithmic scale. 



units at every bias point. Increasing the amplitude h ac 
shifts the smooth branch towards high voltages. In some 
current ranges we observe locking of the I-V curves to 
vertical Shapiro steps at constant voltage. Figures ^a) 
and^Jb) correspond to two different damping constants 
ao = 0.05 and ao = 0.02, respectively. 

The dependence of the voltage 9 7 taken at the con- 
stant dc bias levels 7 = 0.072 and 7 = 0.088 on the 
applied ac-drive amplitude h ac is presented in Fig.^Jc). 
The best linear fits to the power law 0/ cx /if c are with 
the indices /3 = 1.27±0.08 and /3 = 0.85±0.08 for smaller 
and larger current values, correspondingly. 

Thus, we find in good accord with the experimental re- 
sults that the characteristic voltage values of the low volt- 
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age branch depend on the microwave power as V oc P 1 / 2 
(see Fig. 01°) )■ By analogy with the flux-flow behavior, 
we argue that the characteristic voltage V is proportional 
to the average number of fluxons in the junction, which 
in its turn is roughly proportional to the microwave field 
intensity ~ P 1 / 2 . We should take here into account that 
the microwave current generates magnetic field at the 
junction boundaries. At every time moment, the gener- 
ated magnetic field components at the two long junction 
boundaries are equal in magnitude but opposite in po- 
larity. If the field amplitude is larger than critical, there 
will be vortices injected at one junction boundary and an- 
tivortices injected at the other boundary. At microwave 
frequencies / -C f p the magnetic field generated at the 
junction boundaries varies relatively slow in comparison 
with the vortex injection rate. Thus, at each junction 
boundary, there will be a bunch of vortices injected every 
half-period of microwave field and a bunch of antivortices 
injected at every other half-period. During half-period 
of the microwave oscillation, vortices on one side of the 
junction and antivortices on its other side are injected si- 
multaneously and propagate through the junction under 
the influence of the dc bias current. Vortices and an- 
tivortices will eventually annihilate in the middle of the 
junction or, as the dissipation is rather low, may propa- 
gate through one another and annihilate at the opposite 
junction boundary. 
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FIG. 5: Calculated voltage spectra at two bias points of the 
current- voltage curves shown in Fig.f4] (a) point A; (b) point 
B. The numbers near the peaks indicate the corresponding 
harmonics of the ac drive frequency. 

We found that voltage oscillations in the junction 
are very irregular at nearly any point of smooth I-V 
branches. We calculated the Fourier power spectra of 



the voltage oscillations (pit) at the center of the junc- 
tion using FFT. Figure El illustrates such spectra for 
two points A and B of the current-voltage characteris- 
tics shown in Fig. Ufa) . Point A displays a typical spec- 
trum for the states along the smooth branches (see Fig- 
ure EJa)). Voltage oscillations are strongly irregular, the 
oscillation spectrum is very broad and resembles a chaotic 
state. Point B has been chosen at a Shapiro step, which 
is expected to be synchronized with the frequency of ac 
drive / = ui/(2ir) w 0.0318. In the corresponding spec- 
trum presented in Fig. Efb) we clearly see strong peaks 
at frequencies corresponding to various harmonics of the 
driving frequency. We suppose that the prevailing even 
harmonics are related to symmetric injection of vortices 
and antivortices into the junction from two boundaries. 

Next we turn to a numerical simulation of our LSM 
experiment where the laser beam produces a local heating 
of the junction around the beam position at x = xq. 
Thereby the dissipative term a is locally increased. We 
model the laser beam perturbation by a dissipative spot 
placed at the point x = xq and described by the following 
function: 



a(x) = ozq [1 + s6(x — xq)] . 



(4) 



In simulations, we approximate the (5— function by the 
expression [2(j 



eS(x — Xq) ps K 



1 - tanh' 



2(x - x ) 



V 



(5) 



where e — KT). In order to model the laser beam scan- 
ning procedure we calculate the time averaged voltage 
across the Josephson junction as a function of x . This 
algorithm accounts for the experimental situation since 
the period of the fluxon oscillations is much shorter than 
the characteristic times of scanning and power modula- 
tion of the laser beam. Similarly to the experiment, the 
LSM patterns are presented in form of the dependence 
of the inverted voltage change — AQ(xq) at a fixed bias 
point on the I-V curve. In simulations we have chosen 
the strength of the laser spot given by k = 2 77 = 1 . 

Figure shows the calculated voltage response 
—AQ(xo) obtained for the bias point C marked in Fig.f4] 
The scan shown by solid circles has been simulated at 
h = 0. Two types of triangles show simulation results 
for positive and negative dc magnetic field. These plots 
should be compared to experimentally measured ones 
presented in Fig. Elb). One can see that simulations 
qualitatively match the experiment. In both cases, the 
response increases with the magnitude of magnetic field 
and becomes asymmetric. Inverting the field reverses the 
response with respect to the center of the junction. The 
magnetic field makes it favorable for vortices (or antivor- 
tices for negative fields) to enter from one boundary of 
the junction. The response is the largest near the bound- 
ary |l6l Il7| , i.e. in the region where vortices enter the 
junction. 
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FIG. 6: Simulated LSM response of the junction in a state 
corresponding to the marked bias point C in Fig.fJIa), at zero 
(solid circles) and non-zero dc magnetic field (solid triangles 
for positive fields and open triangles for negative fields) . 



IV. DISCUSSION AND CONCLUSIONS 

Our numerical studies show that applying of a mi- 
crowave radiation of large amplitude induces a vortex 
flow in long Josephson junctions. As the frequency of ac 
drive is small, i. e. / <C f p , the dynamics of vortices is 
extremely irregular. It is appropriate to note here that 
chaotic dynamics in long junctions driven by ac field has 
been observed and discussed before @- However, the 
novelty of our results is in the particular appearance of 
the chaotic states in I-V curves, which we relate to exper- 
imental observations. Indeed, the observed low-voltage 
branches are similar to the well-known displaced linear 
branch, which is associated with chaotic dynamics in the 
absence of microwaves 7|- In both cases, the smooth 
branches shift towards high voltages with increasing the 
amplitude of the field (either ac or dc) . This feature be- 
comes even more pronounced, displaying the downward 
curvature (see Figs. EUa) andEJa)) as the chaos in vor- 
tex dynamics is reduced. Note here that the increase of 
temperature (and thus increase of dissipation) suppresses 
the chaos. 

The presence of strong irregularities in the vortex mo- 
tion allows us to draw an analogy with the properties 
of extremely small Josephson junctions subject to mi- 
crowave radiation in the presence of large thermal fluc- 
tuations 0. Indeed, in this case smooth low voltage 
branches in the I-V curve have been observed, and these 
branches also shift towards larger voltages proportionally 
to VP. This effect has been explained as an incoher- 
ent multi-photon absorption by the Josephson phase de- 
scribed by the diffusive motion induced by thermal noise. 
In this case we can write the lumped Josephson phase as 

<p(t) = Vt + ip(t)--cosujt + tp 1 (t) , (6) 
a 

where rj is the amplitude of ac drive, the random function 
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FIG. 7: The power dependence of I-V curves calculated using 
Eq. Q for two values of parameter 8: (a) S = 0.0001; (b) 
S = 0.0002. The power levels correspond to P = 2, 5, 7 dBm. 
In both cases the linear slope corresponding to quasiparticle 
resistance R — 50 kQ (a) and R = 140 kQ, (b), is added. 



ijj(t) determines the Josephson phase diffusion, and (fii(t) 
has to be found self-consistently ^(}. As the frequency 
lu is small, we have obtained the low voltage branch of 
I-V curve in the following analytical form: 

I s (V) = I c ^[f+(V-V s )f + (V+V s )-f-(V-V s )f.(V+V s )} , 

(7) 



where the functions f±{x) = V5^EIEe!e£ 
maximum voltage of the branch V s = 



Here, the 
P, and 



2 +8 2 ' 

rj/a tx 

5 is the parameter characterizing the strength of fluc- 
tuations. The I-V curves are presented in Fig. 7 for 
different values of microwave power and two values of S. 
Note here, that the low voltage branch appears in the 
limit of u> -€1 5 as the thermal fluctuations wash out the 
Shapiro steps. As the fluctuations characterized by the 
parameter S increase, the I-V curve gets more linear (see 



8 



Fig. 7(b)). 

We see that all features obtained for microwave in- 
duced irregular vortex flow in long junctions resemble 
ones observed in a small microwave driven Josephson 
junction in the presence of thermal fluctuations. Thus 
we argue that the ansatz © can be qualitatively used in 
the case of a microwave induced vortex flow, where the 
random function ip(t) is due to chaos in vortex motion. 
The chaos in the vortex motion mimics thermal fluctua- 
tions and wash out the I-V curve. According to both our 
experiment and numerical study, the reduction of chaos 
by increasing the damping (i. e. decreasing 5) leads to 
the pseudo-resonant feature on the I-V curve (see Figs. 
13 a) and0{a)). It is interesting to note that an increase 
of temperature (dissipation) in the case of long Joseph- 
son junction leads to smaller chaos-induced fluctuations. 
This is opposite to the case of a small Josephson junc- 
tions where the parameter 8 increases with an increase 
of temperature. 

In summary, we reported here experimental and 
numerical study of microwave-induced effects on the 
current-voltage characteristics of long Josephson junc- 
tions. We observed that low-frequency microwaves act 



on the junction I-V curves similar to the effect of dc 
magnetic field, i.e. induce the flux-flow behavior. Ex- 
periments using a low-temperature laser scanning micro- 
scope unambiguously indicate the motion of Josephson 
vortices driven by microwaves. Numerical simulations 
are in a good agreement with experimental data and indi- 
cate that the vortex flow is strongly irregular. Our results 
provide also an explanation for previously measured I-V 
curves of intrinsic high-T c Josephson junctions driven by 
microwaves 0- 

When this work was completed, we also learned about 
recent numerical data of the Tuebingen group |2l| . Their 
simulations using the stacked junction model with the pa- 
rameters of intrinsic junctions also show chaotic regimes 
in the vortex dynamics driven by microwaves, which is 
consistent with our observations. 
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